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Abstract 

We study the many body quantum evolution of bosonic systems in the mean field limit. The 
dynamics is known to be well approximated by the Hartree equation. So far, the available results 
have the form of a law of large numbers. In this paper we go one step further and we show 
that the fiuctuations around the Hartree evolution satisfy a central limit theorem. Interestingly, 
the variance of the limiting Gaussian distribution is determined by a time-dependent Bogoliubov 
transformation describing the dynamics of initial coherent states in a Fock space representation 
of the system. 

1 Introduction 

A quantum mechanical system of bosons in three dimensions is described by a wave- function ipN & 
Lg(M^^), the subspace of L^(M^^) consisting of functions which are symmetric w.r.t. permutations 
of the particles. Every physical observable of the system is associated with a self-adjoint operator 
A on the Hilbert space L'^(R^^). The expectation of A in the state described by the wave function 
'i/'Ar is given by the L^-inner product {ipj^,AijjN). For example, the position of the j-th particle is 
associated with the multiplication operator A = Xj. The momentum of the j-th particle, on the 
other hand, is associated with the differential operator A = —iVx^- 

The time evolution of the quantum system is governed by the A^-particle Schrodinger equation 

idt'ipN,t = HN'lpN,t 

where 'ipN,t G L^(]R'^^) indicates the wave-function of the system at time t. On the r.h.s., Hj\[ is a 
self-adjoint operator acting on L^(M^^), known as the Hamilton operator, or the Hamiltonian, of 
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the system. It typically has the form 

N 

j=l i<j 

where A G M is a coupling constant, and y is a two-body interaction among the particles. One could 
also introduce a term describing an external potential (for example, a trapping potential Vextix) — ?• oo 
as |x| — >• oo); the analysis below applies also in this case (under minor assumptions on Ve^t, needed 
essentially to guarantee the self-adjointness of H). 

We are interested in the mean-field regime, where particles undergo a large number of very 
weak collisions, so that the sum of the interactions experienced by each particle can effectively be 
approximated by an average, mean-field, potential. The mean-field regime is realized when ^ 1 
(many collisions), |A| <C 1 (weak interactions), so that A^A =: k is of order one (so that the total 
force acting on each particle is of order one). In order to study the mean-field regime, we analyze 
the dynamics generated by the Hamilton operator 

N 

j=l i<j 

in the limit of large N. In particular, we are interested in the evolution of factorized initial data, 
given by 

N 

V'Ar(xi, . . . , Xn) = JJ (p{Xj) 
J=l 

Although factorization is not preserved by the time-evolution, because of its mean-field character, it 
turns out that, for large A^, factorization is approximately (and in an appropriate sense) preserved, 
i.e. 

N 

V'Ar,t(xi, . . . ,3;Ar) ~ (^j(xj) . (1.2) 

i=i 

A simple argument then shows that ipt must evolve according to the effective nonlinear one-particle 
Hartree equation 

idtipt = -^^^^1 + n{V *\^t\'^) ■ (1-3) 

More precisely, it turns out that (|1.2p can be understood as convergence of the reduced density 
matrices associated with ^N,t- We define the density matrix associated with ipN,t as the orthogonal 
projection jN,t = \i^N,t){'^N,t\ onto ipN,t- Then, for k = 1,...,N, the A;-particle reduced density 
matrix associated with ipN,t is defined as the partial trace of 'yN,t over the last (A^ — k) particles. 

In other words, 7^^^ is defined as a non-negative, trace class operator on L^(R^'^) with integral kernel 
given by 

TyVjiV"^!' • • • ' "^k'-i ^i, ■ ■ ■ , Xf,) 

= j dxk+i . . . dxNlN,t{xi, . . . ,Xk,Xk+i, . . . ,xn;x[, . . . ,x'k,Xk+l, . . . ,xn) 

= j dXk+1 ■ ■ ■ dXNlpN,t{xi, . . . ,Xk,Xk+l, . . . ,XN)lpN,tiXl, ■ ■ ■ ,x'k,Xk+l, ■ ■ ■ ,Xn) 

For a A:-particle observable O^'^) (an observable depending non-trivially only on k particles), we find 
{i^N,t, (O(') (S> 1) i^N,t) = Tr (O(^) 1) -fN,t = TrO('^)7S,l 
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It turns out that the language of the reduced density matrices is the correct language to understand 
the approximation (|1.2p . 



Theorem 1.1 (Theorem 1.1 in [2]). Suppose that the potential V satisfy the operator inequality 
V{x)'^ < D{1 — A), for some constant D > 0, and let H]\f be defined as in Let ijjN = y^^^ , 

i^N,t = e~'^^^^il}N . Then, for every /c G N there exist constants C,K > such that 



Tr 



l^Nl^\^t){^r' <C^ (1.4) 



,K\t\ 



for all t £ M and all N large enough. Here ipt is the solution of the Hartree equation il.3\) . with 
initial data ^t=o = f- 

Remark that the operator inequality V'^ < D{1 — A), which means that 



dxV^{x)\ip{x)\'^ < 



Hi 



for every ip £ L^(]R^), is satisfied by potentials with a Coulomb singularity V{x) = ibl/|x| 
Note that (jl.4p implies that, for any A;-particle observable O^'^^ (X" 1^^"'^), we have 



as — )■ oo, with rate of convergence of order N^'^. In other words (jl.4p implies that, to compute 
the expectation of an arbitrary observable depending non-trivially on a fixed number of particles, 
in the limit of large A^, we can replace the full solution of the A^-particle Schrodinger equation by 
products of the solution of the Hartree equation (jl.3p . 

The first result in the direction of Theorem 11.11 was obtained in [30J for bounded interaction 
potentials, studying the evolution of the reduced densities, as described by the so called BBGKY 
hierarchy (this approach does not give an estimate on the rate of convergence). The approach of 
[30] was then extended to potential with singularities in [131 E El Hffl [Ml HI- More recently, it was 
also extended in [101 \TT\ [12] to the ultradilute limit of quantum mechanics, where particles interact 
through strong potentials with scattering length of the order A^~^; in this case, the many body 
dynamics is asymptotically approximated by the Gross-Pitaevskii equation (mathematical questions 
related with the Gross-Pitaevskii limit were also discussed in [251 13]) • A different derivation of the 
Gross-Pitaevskii equation has been proposed in [28] . 

Another approach to the mean- field limit of quantum mechanics, based on ideas from [191116] has 
been obtained in [29] and later extended in |261 [2] . This approach is based on a representation of the 
system on the bosonic Fock-space and on the use of initial coherent states and it has the advantage 
of giving precise bounds on the rate of convergence towards the Hartree dynamics. Theorem ll.il as 
stated above, has been proven in [2] making use of these techniques. The approach of [29] is also the 
basis for the analysis in the present paper; we will therefore describe it in details starting in Section[2l 
Alternative approaches and questions related with the mean-field limit of quantum dynamics have 
been discussed in [9l EH H [fsl E] . 

Because of the bosonic symmetry, physically measurable observables on L^(R'^^) are invariant 
w.r.t. permutations of the A'^ particles. A physically measurable one particle observable has the form 
O = Ylf=i 0^^\ where O^^^ = Kg) ■ ■ ■ O 010 ■ ■ ■ 01 acts non-trivially only on the j-th particle, 
for a self-adjoint operator O on L^(M^). The kinetic energy '}2f=i is an example of such an 

observable. The operator X^^Li XA{xj), measuring the number of particles in a set ^ C is another 
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example. In probabilistic terms, one can think of the observables O^^^ as random variables, whose 
distribution is determined by the wave function ijjN ^ L^(R^^) describing the state of the system (the 
distribution is determined through the spectral decomposition of the observables). If V'at = , the 
are independent identically distributed random variables with a common distribution determined 
by the one-particle wave function ip G L^(M'^) and, as a sum of independent and identically distributed 
random variables, O = XljLi O^'^^ satisfies a law of large numbers and a central limit theorem. If we 
let the system evolve, ipN,t = e~^^^^ip^^ does not factorize. Hence, at time t ^ Q, the variables O*--'-* 
are not independent (they are still identically distributed because of the permutation symmetry of 
fpN,t)- Nevertheless, Theorem 11.11 implies that, asymptotically, iljN,t can still be approximated in a 
certain sense by a factorized wave function. As it turns out, the result of Theorem 1 1 . 1 1 easilv implies 
that, for any time t G M, O = Ylf=i O^'^^ still satisfies the law of large numbers (for a bounded 



self-adjoint operator O on 



In fact, let O = O 



O — {ipt,Oipt), and let 



and 



"<PN,1 



denote probabilities, respectively, expectations w.r.t. the distribution determined by the wave 



function tpN f Then, by Markov's inequality. 



1 ^ ~ 



N 



>e < 



From Theorem I l.H we conclude that 



N 



^Tr7S,'i(0»0) + ^T: 7)^,^0 



(i)n2 



lim sup P^j^^ 

N^oo 



1 ^ ~ 

N ^ 



In the next sections, we show that, for all times t E 
a central limit theorem, in the sense that, as — )• oo 



the one-particle observable O also satisfies 



1 



N 



converges in distribution to a centered Gaussian random variable. The variance of this Gaussian 
variable is different from the one we would observe replacing the solution ■i/'Tv,* of the Schrodinger 
equation by the product ip®^ (the correlations among the particles developed by the A^-particle 
dynamics are sufficiently weak to still have a central limit theorem, but they are sufficiently strong 
to affect the variance of the limiting Gaussian fluctuations). As we will see the limiting variance is 
determined by a so called Bogoliubov transformation which describes the dynamics of fluctuations 
around the mean fleld Hartree dynamics in a Fock space representation of the system, which will be 
deflned in the next section. 

This is not the flrst time a central limit theorem is proven in a quantum setting. Quantum central 
limit theorems stating the convergence of polynomials in non-commuting self-adjoint observables to 
appropriate Gaussian distributions already appeared in the 70's, in the works [UEO]. Since then, these 
results have been extended and applied to questions in equilibrium quantum statistical mechanics 
in [T7j, non-equilibrium quantum statistical mechanics in [21], quantum information theory in |18j . 
and combinatorics in [22j. The paper ^ is closer to the spirit of our manuscript, since it deals with 
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quantum dynamics and with non-independent variables; however, in our case the correlations are 
generated by the time evolution (the particles are uncorrelated at time t = 0), while in [5] they are 
already present in the initial state (the Hamiltonian in this case is quadratic and does not generate 
correlations for factorized initial data). 

2 Fock space representation and fluctuation dynamics 

To state the central limit theorem, we first have to study the fluctuations of the many-body evolution 
around the limiting Hartree dynamics. To this end, it is convenient to introduce a Fock-space 
representation of the system under consideration; this will allow us to study the evolution of so- 
called coherent states, for which the fluctuations can be written in a nice and compact form. 
The bosonic Fock-space over L^(M^) is defined as the direct sum 

T = C®^L'l{M.^'',dxi,...dxn) 

where L^(]R^") denotes the subspace of L^(R'^") consisting of functions symmetric with respect to 
permutation of the n particles. Vectors in the Fock-space are sequences n>o, where 

is an n-particle bosonic wave function. The idea behind the introduction of the 
Fock-space is that we want to study states with non fixed number of particles. Clearly, has the 
structure of a Hilbert space with the inner product 

n>l 

The vector = {1, 0, 0, ... } G is called the vacuum and describes a system with no particles. An 
important operator on is the number of particle operator AA, which is defined by 

The vacuum S7 is an eigenvector of M with eigenvalue zero. More generally, vectors of the form 
{0, . . . , 0, , 0, . . . } , having a fixed number of particles, are eigenvectors of Af (with eigenvalue m) . 
The number of particle operator is an example of second quantization of a self-adjoint operator on 
the one-particle space L^(M'^). In general, if O denotes a self-adjoint operator on L^(M'^), we define 
its second quantization dr{0) as the self-adjoint operator on T defined by 

n 

where O^^^ = 1 (8) • • • O ® • • • (8> 1 acts as O on the j-th particle and as the identity on the other 
(n — 1) particles. With this notation M = dT{l). 

On J^, we define the Hamilton operator TIn by ?^7v{V'^"^}n>i = {'HS^V^"^} 

„>i, with 

n 1 " 

= J2 -A., + ^ E ^(^^ - ■ (2.1) 

j=l i<j 

By definition, T-Ln leaves each n-particle sector J^n (defined as the eigenspace of M associated with 
the eigenvalue n) invariant. Moreover, on the A^-particle sector, T^^r agrees with the mean field 
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Hamiltonian Hi\i defined in (jl.ip . In particular, if we consider the Fock-space evolution of an initial 
vector with exactly particles, we find 

e-'^^'{0, . . . , 0, Vtv, 0, . . . } = {0, . . . , 0, e-^^^ViV, 0, . . . } 

exactly as in Section [TJ The advantage of working in the Fock space is that we have more freedom in 
the choice of the initial data. We will use this freedom by considering a class of initial data, known 
as coherent states, with non-fixed number of particles. 

It is very useful to introduce, on the Fock space J-', creation and annihilation operators. For 
/ G L^(M'^), we define the creation operator a*{f) and the annihilation operator a(/) by setting 

1 " 

{Xl, . . . , Xn) = ^ ^ /(Xj)V'("~^)(xi, . . . , Xj-l,Xj+l, ...,Xn) 

^=1 (2.2) 

(a(/)V')(") {xi,...,Xn) = V^^ I dx7MV'^"+'^(a:,a;i,...,x„). 



The operators a*{f) and a{f) are densily defined and closed. It is easy to check that, as the notation 
suggests, a*{f) is the adjoint of a{f). Creation and annihilation operators satisfy the canonical 
commutation relations 

[a{f),a*{g)] = {f,g) and [a{f),a{g)] = [a* {f) , a* {g)] = Q (2.3) 

for any f,g £ L^(M^) (here {f,g) denotes the L^-inner product). Physically, the operator a*{f) 
creates a particle with wave function /, while a(/) annihilates it. As a consequence, a state with N 
particles all in the one-particle state cp, can be written as 

{o.^^.o.^«",o,...} = <5^n. 

It is also useful to introduce operator- valued distributions ax,a* defined so that 

a{f) = j dxf{x)ax, and a*(/) = jdxf{x)al. (2.4) 

With this notation, a%ax gives the density of particles close to x S M^. The number of particles 
operator can formally be written as 

M = j dx a* ax ■ 

More generally, the second quantization of a self-adjoint one-particle operator O on L^(M^) with 
integral kernel 0{x,y) can be expressed as 

dT{0) = J dxdy 0{x,y)al.ay . 

Similarly, the Hamilton operator H.^ can be formally expressed in terms of operator- valued distri- 
butions as 

'Hn = j dxVxal VxO-x + ^ y ^^^y ~ y)ala*yayax . (2.5) 

The fact that every term in the Hamiltonian contains the same number of creation and annihilation 
operators means that T-Ln commutes with the number of particles or, equivalently, that the number 
of particles is preserved by the time-evolution enerated hy Hn- 
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For later use, we observe that the creation and annihilation operators are not bounded; however, 
they can be bounded by the square root of the number of particles operator, in the sense that 

\Hm\ < 11/11 IIAAi/^V^II and ||a*(/)V^|| < ||/|| ||(AA + 1)1/^11 (2-6) 

for every il) ^ F , f ^ L^(M^) (here ||/|| indicates the L^-norm of /). 

As mentioned above, we are going to study the evolution of initial coherent states. For arbitrary 
if G we define the Weyl operator 

Ty((^) = e"*(^)-"('^). (2.7) 

The coherent state with wave function (p is then defined as W{ip)Q.. The Weyl operator W{ip) is a 
unitary operator; therefore the coherent state W{'^)^ always has norm one. Moreover, since 

the coherent state W{(p)^ does not have a fixed number of particles. One can nevertheless compute 
the expectation of the number of particles in the state W{ip)rL] it is given by 

{w{^)n,Mw{ip)n) = . 

More precisely, it turns out that the number of particle in the coherent state W{(p)VL is a Poisson 
random variable with expectation and variance ||</?|p. The main reason why coherent states have 
nice algebraic properties (which will be used later on in the analysis of their evolution) is the fact 
that they are eigenvectors of all annihilation operators. Indeed 

a{f)W{ip)^ = {f,ip)W{ip)^l 

for every f,(p£ L^(M^). This is a simple consequence of the fact that Weyl operators generate shifts 
of creation and annihilation operators, in the sense that 

W*{ip)a{f)Wi^) = aif) + if, if), and W*{ip)a*{f)W{ip) = a*{f) + {if, /) . (2.8) 

Instead of studying the time evolution of initially factorized states, of the form = a*{(p)^Q,/^/W., 
it is convenient to analyze the dynamics of coherent states of the form W{\nV(p)i}, for (f E L'^(M^) 
with 1 1 (/J 1 1 = 1 (at the end, it is possible to express factorized states as linear combinations of coherent 
states, and therefore to establish properties of their evolution as well). The expected number of 
particles in the state W{\/T^ip)Q, is equal to A^. For this reason, we may expect the evolution of 
W{^/N(p)^} to have the same mean-field properties as the one of ipN. In other words, we may expect 
that the evolution of W{\nV(p)i} is still approximated by the coherent state W{^/Nipt)0,, where cpt 
is the solution of the Hartree equation ()1.3p . with initial data (pt=o = f- To study the fluctuation 
around the mean-field approximation H^(\/iV(/7t)f], we introduce the two-parameter group of unitary 
transformations 

U{t; s) = e-'^^*''^ T4^*(\/iV^t)e-*^^(*-^%(\/iV^,) (2.9) 
for all t, s G M, with the phase factor 

^{t'^s) = — J J dx{V *\(fr\^){x)\(pr{x)\'^. 
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It turns out that V({s; s) = 1 and that U{t; s) is strongly difFerentiable on the domain P(/C +M) C -F 
(at least under the assumptions on (p and V given in Theorem II. ip . Here /C is the kinetic energy 
operator 

IC = dT{-A) = j dxV^a*V^a^ (2.10) 

The derivative of U{t; s) is then given by 

idtU{t]s) = C{t)U{t;s), (2.11) 

where we introduced the generator 

C{t)= / Va;a*Va;aa; + / dx {V * \ipt?){x)alax + / dxdyV{x-y)(pt{x)lpt{y)alay 



(2.12) 



+ dxdyV{x-y) [Lpt{x)ipt{y)ala*y + Lpt{x)ipt{y)axay) 
+ j dxdyV{x - y) a* (a*y?t(y) + ayTpt{y)) 



The phase factor Lo{t]s) in (12. 9p really does not play an important role (it cancels from expressions 
of the form U*(t; s)AU(t; s), for any operators A on J^). To understand the effect of the unitary 
evolution lA{t;s), we observe that, for a one-particle operator O on L^(M^) with integral kernel 
0{x,y), we have 

= ^ j dxdyO{x,y) (^e-''^^^W{^/Nip)n, a^ay e-''^^*W{VN^)n 

_ 1 

~ N 



dxdyO{x, y) {U{t; 0)n, (a* + VNi^t{x)){ay + V N ipt{y)) U {t; 0)n 



Coo{t)= j dxV^alVxa^ + / (F * p)(x)a*a^ + / dxdyV{x-y)iptix)(pf-{y)al.a 
+ \ I dxdyV{x-y) {iptix)ipt{y)ala*y +Tpi{x)i^t{y)axay) 



Hence, U{t]s) describes the evolution of the fluctuations around the mean-field dynamics. As — ?• 
oo, the last two terms on the r.h.s. of (I2.12P vanish. We can therefore expect the fluctuation dynamics 
hl{t] s) to converge in this limit towards a limiting evolution with time-dependent generator 



(2.13) 



The existence of such a limiting dynamics was established in [16]; for completeness, we reproduce 
the relevant statement in the next proposition. 

Proposition 2.1 (Prop. 4.1 in [16]). Suppose V £ L°°{R^) + L'^{M.^), and that the map t (pt 
is in C(]R, L-^(M^) n L^(M^)) (both conditions are easily checked under the assumptions of Theorem 
\3.1\ below). Then there exists a unique two-parameter group of unitary transformations Uoo{t', s) with 
i^oo{s; s) = 1 for all s S M, and such that Uoo{t] s) is strongly differentiable on the domain 'D{K, +M) 
(7C is the kinetic energy operator from i2.10\) ) with 

idtUooit; s) = Coo{t)Uoo{t; s) . (2.14) 
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It turns out that the Hmiting dynamics Uoo{t; s) can be described as a so called Bogoliubov trans- 
formation. This is a consequence of the fact that the time-dependent generator Coo{t) is quadratic 
in the creation and annihilation operators. For f,gG L^(]R^), we introduce the notation 

Aif,g)=aif) + a*ig) 

We observe that 

{A{f,g)r =A{gJ) = A(^(^^j o)^-^'^^ ^^'^^^ 

where J : L^(M^) — )• L^(]R^) is the antilinear operator defined by Jf = f. From p.3p . we find that 
the operators A(f, g) satisfy the commutation relations 

[A{h,g,),A*{h,g2)] = {{fugi),S{f2,g2))L^^L^ with 5 = J \ ^ (2.16) 

Theorem 2.2. Assume V"^ < D{1 — A) and if E i/^(M^). Then, for every t,s E M there exists a 
bounded linear map 

e{t; s) : l2(r3) e l2(r3) ^ l2(m^) e ^^(m^) 

such that 

U*^{t; s)A{f,g)U^{t; s) = A{e{t; s){f, g)) 
for all f,g£ L^(M'^). The map @{t; s) is so that, for every t, s E M, 

o)=(j 0^e(,;,) (2.17) 

and 



S = @{t; s)* S @{t; s) (2.18) 

Moreover, Q{t; s) is such that 

@{t;s){ipt,lpt) = {y^s,^s)- (2.19) 

A map Q{t;s) with the properties {2.11^ and i2.18\) is known as a (hosonic) Bogoliubov transfor- 
mation. Note that \2.11^ and 112. 18\) simply state that i2.15]) and, respectively, 112. 1 6\) are preserved 
under Q{t]s). Like every Bogoliubov transformation, 0(t;s) can be written as 



V{t;s) JU{t;s)J 

for bounded linear maps U{t;.s),V{t;.s) : L^{R^) L^{R^) satisfying 

U*{t; s)U{t; s) - V*{t; s)V{t] s) = 1 and U*{t; s)JV{t; s)J = V*{t; s)JU{t; s)J. 

The maps U (t; s) , V (t; s) are not only bounded from L^(M'^) into L^(M'^). Instead they even extend 
to bounded maps from H^iW^) into H^(W^) (and @{t;s) extend therefore to a bounded maps from 
H\M.^)eH^{M.^) into itself). 

Remark. The Bogoliubov transformation Q{t; s) is a two-parameter group of symplectic trans- 
formations. Formally, they satisfy the equation 

idte{t] s) = e{t; s)A{t) (2.21) 
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with generator 

V Bt -JDtJ 
where we defined the hnear operators 

Dtf = -A/ + (V * \^t\^)f + {V * ^J)ipt 
Btf = {V*Tptf)Tpt 

acting on L'^{R^). We observe that = Dt, = JBtJ. and therefore Al = SAtS {A{t) is self- 
adjoint w.r.t. the degenerate product defined by S; hence Q{t; s) is unitary w.r.t. this inner product, 
which means that @{t;s) is symplectic) . Note also that the condition V'^ < D{1 — A) implies that 
Bt is a Hilbert-Schmidt operator for all t G M. This is an important fact to make the Bogoliubov 
generated by A{t) implementable. Since the existence of a generator A{t) for @{t;s) is not needed 
for our analysis, we do not try to make ()2.2ip precise. 

We will prove Theorem 12.21 in Section HI where we establish several properties of the limiting 
evolution Uoo)t; s). 



3 The Central Limit Theorem 



We are now ready to state our main result. 

Theorem 3.1. Suppose V is a measurable function such that V'^{x) < D{1 — A), for some D > 0. 



For ip £ 



with \\ip\\ = 1, let •i/'AT = and ^N,t 



">pN- Let O he a hounded self adjoint 



operator on L^(M^) with ||VO(l — A) '^^'^\\ < oo. For arbitrary fixed t £ 



let 



N N 
i=i i=i 



(3.1) 



where O^^'^ = l^-?"^) ® O ® i(^-i) denotes the operator O acting on the j-th particle. Then, as 
N — )• oo, the random variable Ot (whose law is determined by the N -particle wave function ipN,t) 
converges in distribution to a centered Gaussian random variable with variance 



{Bit- 0) {Oift, JOift) , G(t; 0) {Oipt, JO^t)) 



e(t;0)(O(^t,JO(^t),-^ {^,^) 



(3.2) 



\\u{t-d)o^t + jv{t-^)o^t 



\{^,U{t-Q)0^t + JV{t-Q)0^t)\^ >Q 



where the bounded linear operators U{t;Q),V{t;^) : 
L2(m3) _^ l2(^m3) © l2(^m3) are defined in Theorem\2M 



and e(t;0) : L'^( 

Notice that, for t = 0, since 0(0; 0) = 1 (and thus f/(0;0) = 1, y(0;0) = 0), the variance is given 



by 



which is consistent with the fact that, at time t = 0, V'Af,t=o = 'p'^^ and the variables O^^^ are there- 
fore independent identically distributed random variables, with a common distribution determined 
(through the spectral theorem) by if. 

Instead of the variable Ot, we could also consider the variable 

N N 
.7 = 1 .7 = 1 
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which is centered for ah N (while Ot is only centered in the limit — )• oo). Observe that 



0[ = 0t + Sn, with 6n = ViVTrO {\ipt){vt\ - tS 

It follows from Theorem 11.11 that, for every fixed time, \6n\ ^ N^^^"^. For this reason, like Ot, also 
the random variable O't converges in distribution, as — t- oo, to a Gaussian random variable with 
variance erf as given in (j3.2p . 

We also remark that the assumptions (p S H'^{M^) and ||V0(1 — A)""*^/^!! < oo are needed to 
control the singularity of the potential V. If instead of V'^{x) < D{1 — A) we just assumed that 
the operator norm of V is bounded (i.e. that the L°° norm of the function V is finite), the result 
would hold for arbitrary if £ H^{W^) and for arbitrary compact operator O and its proof would be 
substantially simpler (the analytic part of the proof would be easier, the combinatorial part would 
be unchanged). Since however we believe the Coulomb potential to be physically interesting, we 
preferred to avoid the assumption of bounded potential. 

For example. Theorem 13.11 applies when O = xa{x), where XA is a smoothed out version of the 
characteristic function of a set ^4 C M'^ (one needs the smoothed out version, because of the condition 
\\Vxa{1 - A)"i II < oo). In this case, the operator J2i=i XA{xj) measures the number of particles in 
the set A. The law of large numbers implies that X^^Li XA{xj) converges almost surely to 

N{^t,XA{x)^t) = N J dxxAix)\Mx)\^ (3.3) 

as — )• oo. Theorem 13.11 implies that the typical fluctuations around (j3.3p are of the order A^^/^ 
and that, on this scale, they approach a Gaussian law. 

In order to prove Theorem 13.11 we show that the moments of the random variable Ot, with 
distribution determined by the wave function ipN,t converge, as A^ — )■ oo, to the moments of a 
centered Gaussian random variable with variance erf. 

Lemma 3.2. Under the assumptions of Theorem \3.1\ we have, for any fixed £ £ N, 

I ,7 \ I if i is odd 



4 Key bounds on fluctuation dynamics 

To show Lemma 13.21 we need some crucial estimate on the fluctuation dynamics U{t;s) defined in 
(j2.9p and on its formal limit Uoo{t; s) defined in (j2.14p . In particular, we need to control the growth 
of the expectation of the number of particles operator M, of its high powers , and of the kinetic 
energy operator 

fC = J dxV xCL*V xCLx 

We assume that the potential V entering the definition of U{t;s) and of Uoc{t;s) satisfies the 
operator inequality 

y2(x) < D{1 - A) (4.1) 

for some D > 0. Moreover, ipt is the solution of the Hartree equation ()1.3p corresponding to an 
initial data ip G if^(M^), with \\ip\\ = 1. Note here that the bounds on the growth of the number 
of particles only depend on the norm of the initial data and therefore can be extended to all 
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Lf G /7^(M^). The bounds on the growth of the kinetic energy, on the other hand, depends on \\'pt\ 
From the Hartree equation (|1.3p . we observe that 



W^tW < IIA^Jtll + \\{V* \ipt\'^)^t\\ < C\\ipt\\H2 < C\Mh2 (4.2) 

uniformly in t (the constant C depends only on the constant D in (|4.ip ). Here, we used the propa- 
gation of regularity for solutions of the Hartree equation (jl.3p : for any s > 1, there exists a constant 
C such that 

supWcptWH" < CWipWH"- (4.3) 

teM 

We will use this fact in several occasions. 

In the next proposition, we control the growth of powers of the number of particles, of the kinetic 
energy and of its square with respect to the limiting dynamics Uooit; s) defined in (|2.14p . 

Proposition 4.1. For every j £ N, there exist constants C,K > (depending on D, \\ip\\iji, and j) 
such that 

{U^{t;s)i^,M^U^{t;s)^) <Ce''^'~'^{i^,{M + iyi>) (4.4) 

for all ip £ T, and all t, s G M. Moreover, there are also constants C,K > (depending on D, 
\\(p\\Hi) andC',K' (depending on D, and \\ip\\ff2) such that 

{Uoo{t;s)iljXUoo{t;s)^l,)<Ce^\'"\{ilj,{lC+M + l)'^) (4.5) 

and 

{Uoo{t;s)i;,}C^Uoo{t;s)i;) <C'e^'\'-'\{iP,{)C^ +Af^ + l)i;) (4.6) 
for all ij: £ , and all t, s G M. 

Proof. Eq. (j4.4p is taken from [2] [Prop. 5.1]. To prove (j4.5p we combine Lemma 14.41 below, where 
we prove the bound K, < Coo{t) + C{N +1) with the estimate 

\{u^{t- s)iij, Coo{t)Uoo{t; s)V')l < Ce^i*-^i(^, {Coc{s) +M+ 

on the growth of the expectation of Coo{t), proven in [2] [Lemma 6.2]. To prove ()4.6p . we write 

}C = Coo{t) - Ai^t - A2,t - As^t 

where 



Ai^t= J dx{V*\(ft\ ){x)alari; 

A2,t = J dxdyV{x - y)Lpt{x)lpi{y)alay 



and A:it = Bt + B^, with 



Hence, we have 



B= dxdyV{x - y)(ft{x)ft{y)ala* 



JC^<C^{tf + Al, + Al, + Al, (4.7) 
Observe that, because of the assumption V'^{x) < D{1 — A), and because of (14. 3D . 

\Ai^t\ < sup \V *\ipt\'^\ M < CM (4.8) 
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for a constant C > depending only on D and on Similarly 



{'^,M,t^)< J dxdy \V {x - y)\\<pt{x)\\ipt{y)\ \\a^ij\\ Wttyi'W 

< / dxdy\V{x - y)\'^\ipt{x)\'^\\ayijf + / dxdy \(ptiy)\'^\\a: 



(4.9) 



< 



iDyt\\Hi + i){^p,M^p). 



Combining the last bound with the corresponding lower bound (and with ()4.3p ). we find \A2^t\ ^ CM. 
Since Ai t and A2 t commute with AA, it follows that 



Al,t < CM^ and Al t < CM^ . 
We still have to consider the term ^3,4. To this end, we notice that 

Alt<BtBl + B:Bt 

where 



(4.10) 



(4.11) 



BtBl = J dxdydzdwV{x - y)V{z - w)(ft{x)<ft{y)<ftiz)^Pti'w)ala*aza^ 
< / dxdydzdwV'^{x - y)\(ft{x)\'^\(pt{y)\'^ a*,a*^ay,az 



(4.12) 



and 



B^Bt = BtB 



f* + 4 J dxdydzV{x - y)V{y - z)\(pt{y)\'^(pt{x)(pt{z)alaz 

+ 2 J dxdyVHx-y)\ipt{x)\^\My)\^ 
< C{Af+lf 



(4.13) 



for a constant C depending only on D and on the norm H^jUji^i. From (|4.7p . (j4.10p . (j4.12p . (j4.13p . 
we conclude that 

<Cl^{t) + C{Af + lf (4.14) 



Similarly, one can also write Coo{t) = IC + Ai^t + ^2,t + ^3,t; and prove that 

Cl,it)<IC^ + C{M + lf 



(4.15) 



Next, we need to control the growth of the expectation of £^(t) w.r.t. the dynamics Uoo{t]s). 
To this end, we remark that 

^ {U^{t; s)ij, Cl^{t)Uoo{t] s)V> = (U^{t] s)V, {^C^{t)C^{t) + Coo{t)Coo{t)) U^{t- s)V 

Hence 



-{U^{t-s)ij,Cl,{t)U^{t-s)^) 



<aA^{t-,s)^,Li,{t)U^(t-,s)^) ' {Uoo{t;s)i;,Coo{tyUoo{t;s)'^ 
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where 

Here we introduced the notation 



^i,t = j dx{V* {iftipt + iftipt)) (x) ala^ 

A2,t= / dxdyV{x - y) {ipt{x)Tp^{y) + ipt{x)Tp^{y)) a^ay 



and ^3,4 = Bt + with 

It follows that 
We observe that 



= 2 y dxdyV{x - y)(ft{x)ft{y) a^a* 

J^lit)<Al, + Al, + Al, (4.16) 



sup\V*{(ptft + fm)\<'^siip I dy\V{x - y)\\(pt{y)\\^t{y)\ 

<2\\ipt\\2 sup( / dyV\x-y)\ipt{y)f 



(4.17) 



for a constant C depending only on D and Hv^Hi^a (here we used (14. 2p ). Eq. (I4.17[) implies, similarly 
to (j4.8p . that Af^ < CAf^ (this time, however, with a constant C depending on the H^-norm of 
the initial data if). Analogously to (14. 9p . we find A^^ < CA/"^ (again with a constant depending on 
||(/9||j:^2). Finally, to bound the contribution from A^^, we observe that 

B B* = / dxdydzdwV{x — y)V{z — w)(f't{x)'Pt{y)^t{^)^t{'^)^*xO'v^z0.w 



<2j dxdydzdwV{x - y) \(pt{x)\ \'Pt{y)\ alal,auiaz 

<\Wt\?H^ \mlN^<CN^ 

for a constant C depending on ||(/7||j:^2. Similarly, we also obtain B*B < C{J\f + 1)^ (the +1 factor 
takes into account terms arising from various commutators; this computation is similar to (|4.13|) ). 
We conclude that 

Cl{t)<CiM+lf 

and thus, with (03]), 
d 



j^{U^{t-s)tlj,Cl,{t)U^{t-s)il^) 



1/2 



Gronwall's inequality implies that 

{Uooit■,s)^l^,CUtpoo{t■,s)i^) < Ce^'*-^l (V^, (£^(0) + + l) ^) 

for appropriate constants C,K > Q. From ()4.14p and ()4.15p . we conclude that 

(Z^oo(t;s)V',/c2z^oo(t;s)V'> < Ce^l*-^l (V-, (£^(0) + i) ^) 

<Ce^l*-^l(V,(^'+AA2 + l) ^) 

□ 
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Proposition 14.11 can be used to prove properties of the Bogoliubov transformation @{t; s) defined 
in (j2.20p . In particular, we can now show Theorem 12.21 



Proof of Theorem It is proven in [2] [Lemma 8.1] that 

PkW^{t;s)a*{f)Uoo{t;s)n = 
PkU*^{t;s)aif)U^{t;s)n = 

for all A; 7^ 1 (here Pk denotes the orthogonal projection on the k-th sector J^^ of the Fock space). In 
other words, the vector l/(^{t; s)a* {f)Uoo{t; s)^l lives in the one-particle sector, for all s G M. This 
implies that there are linear operators U{t; s), V{t; s) : L^(]R^) — J- L^(M'^) with 

U*^{t; s)a*{f)Uoo{t; s)n =a*{U{t; s)f)n 
U*^{t; s)a{f)Uo,{t; s)n =a*{JV{t; s)f)n 

Recall here that Jf = f. The operators U{t; s) and V{t; s) are bounded because, by Proposition 14. H 

\\U{t;s)f\\ = \\a*iUit;s)fM = \\a* {f)U^it; s)n\\ < \\f\\ \m + l)'/^Uooit; sM < Ce^l*-^l ||/|| 
and, similarly, 

\\V{t;s)f\\ = \\a*{JV{t;s)f)n\\ = \\a{f)Uoo{t; s)n\\ < \\f\\ WM'/'^Uooit; s)n\\ < Ce^'*"^! ||/|| . 

Defining the linear map e(t; s) : L'^{R^) L^{R^) L^{R^) L^iR-^) as in ^^2Q^), we obtain 

U*^{t; s)A{f,g)Uoo{t; s)n = A{e{t; s){f, g))n (4.18) 

for all f,g^ L^(]R^). Note that the boundedness of &{t; s) follows from the boundedness of U{t; s) 
and V{t;s). We notice also that the boundedness of U{t;s) and V{t;s) from i/^(M^) into H^(R^) 
follows similarly as above, but using the estimate ()4.6|) for the growth of the kinetic energy. 

Next, we define, for fixed ip G P(/C + M) (this assumption guarantees that Uoo{t;s) can be 
differentiated; see Proposition 12. ip . s G M, (7 € Lp'iR^), and for any bounded operator M on with 
MV{1C + M) CV{}C + M), 



F{t):=Y^ sup ■ 



[[W^it;s)Jif)U^it;s),a}'{g)] ,M] 



where a", a" are either annihilation or creation operators. We observe that, since e *'^*a^(/)e 



, we also have 
[[U^{t;s)e-''^('-'^J{f)e"^('-'^U^{t;s),a^{g)] ,M] ^p\\ 



att(e'*^/), and since ||e**^/|| = 
F{t):=Y, sup 

« /eL2(R3) 



A simple computation shows that F{s) = 0. Moreover, we find 



d 
di 



W^it-s) [(£oo(t)-/C),e-^'=(*-^)a(/)e^^(*-^)l Uo^it; s), a\g) 



W^it-s) r(£oo(t)-/C),a(e-^^(*-^)/)l Uo.it; s) , a\g) 
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Using the canonical commutation relations, and the notation f{t) = e **^/, it is simple to check 
that 

[£oo(t) - /C, a{f{t - s))] = aiiV * \^t\^)f{t - s) + {V * f{t - s)^t)^t) + a*{2{V * J{t - s)ipt)ipt) • 
Hence, we conclude that 



< / dT 

ft 



W^{t- s)a [{V * \iprf)f(T -s) + {V* f{T - s)-^^)^r) U^{t- s),a\g) 



.M 



dT 



U*^{t; s)a* (2(F * /(r - s)ipr)iPr) U^{t- s),a\g) 



M 



<C j'dr {\\{V * Wr\')f{T -S) + {V* f{T - S)ip,)cp^\\ + \\{V* f{T - s)ipr)Vr\\) F{t) 

Notice that, under the assumption V"^ < D{1 — A), we find 

\\{V * \^r\')fir -s)\\< ||/(r - s)\\ sup / dyV{x - y)\My)\' < < C||/|| 

and 

\\{V*fiT-s)lp,)^r\\ < WVrW sup / dy V {x - y)\f (t - s , y)\\ipr{y)\ < C\\f\\ yr\\%i < C\\f\\ 

X J 

for a constant C, independent of t and /. Therefore 

\\[[U^{t;s)e-"^('-'^a{f)e"^('-'^U^{t;s),a^{g)] , M] ip\ 



< C / drF(r) 



The term with a*[f) instead of a(/) can be bounded similarly. Therefore, taking the supremum over 
/ e L^(M^), we conclude that 



F{t) < dTF(T) 

Since F[t) > for all t G M, and F[s) = 0, we conclude that F{t) = for all t G M. This proves that 

[[U*^{t-s)A{h,gr)U^{t-s),A{h,g2)],M\=^ (4.19) 

for every fi, f2, 9i, 92 S L^(M^), and every bounded operator M on such that MV{K, + M) C 
V{1C+M). This implies that 

{^l;,[U*^{t;s)A{h,gi)U^{t;s),A{f2,g2)]i^) = {n,[U*^{t; s)A{fi, gip^dt; s), A{f2, g2)]n) (4.20) 

for all ^ G V{K.+M) with HV'H = 1. This follows because, from (lOTIl . 

[[U*^{t-s)A{h,gi)U^{t-s),A{f2,g2)],P^\ = [[U*^{t-s)A{h,gi)Uoo{t-.s),A{f2,g2)],Pn] = ^ 
where P^,Pn denote the orthogonal projections into ^ and, respectively, Q. Therefore, 

{^^,[U*^{t;s)A{fl,gl)U^{t■s),A{f2,g2)]n) = {i;,[U*^{t;s)A{^,gl^^^ 

= {iP,[W^{t-s)A{h,g^)U^{t-s),A{f2,g2)m 
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and, similarly, 

{i;,[W^{t;s)A{fi,g,)Uoo{t;s),A{f2,g2)] 5^) = 0',^) {n,[U*^{t; s)A{h, gi)Uoo{t; s), Aih, 92)] 5^) 

If 7^ 0, this implies immediately (|4.20p . If = 0, (|4.20p follows repeating the argmnent 

with ^ = 2-1/2 (V^ _^ Combining (^jTM with (^TIM . we find 

{^^,[U*^{t;s)A{h,g,)U^{t;s),A{f2,g2)]^|J) = {n,[A{Q{t■s){h,g,)),A{f^^ 

= (e(t;s)(/i,5(i),5(/2,5'2))L2g5i2 

where S is defined as in ()2.16p . Hence 

[U*^{t;s)A{fi,gi)U^{t;s)-A{e{t;s){fi,gi)),A{f2,g2)] = (4.21) 

for all /i,5i,/2,52 e ^'(K^). Let R = U^{t; s)A{h, gi)Uoo{t; s) - A{e{t; s){h, gi)). We already 
proved in (j4.18p that RQ = 0, and in (j4.2ip that R commutes with any creation and annihilation 
operator. Since states of the form a* (hi) . . . a*(/„)r2 form a basis for J^, this implies that R = 0, and 
therefore that 

U*^{t;s)A{f,g)U^{t;s) = A{e{t;s){f,g)) 

for all f,g e L'^{R^). From 

iAiQ{t;s)if,g))r = (KUt; s)A{f, gp^^it; s))* 
= W^{t-s)A*{f,g)UUt;s) 

= W^{t- s)A{Jg, Jf)U^{t- s) = A{Q{t- s){Jg, J f)) 

we deduce (12.17p . The property (12.18j) . on the other hand, follows because 

[A{Q{t;s){h,gi)),A{Q{t-s){f2,g2))] = [Ul,{t-s)A{h,gi)U^{^^^ 

= uut; s) [^(/i, gi), ^(/2, 52)] Uoo{t; s) 

= ((/l,5'l),5'(/2,fi2))i2g3i2 

Finally, to prove (j2.19p . we observe that 

A{Q{t-s){^uTpt))=W^{t-s)A{^u^t)Uoo{t;s) (4.22) 

We have 

i-^l^^it;s)A{ipt,Tp^)Uooit;s) = -U^{t;s)[Coo{t),A{ipt,Tpt)]Uoo{t;s) 



+ U^{t; s)A{iift, -i'^t)Uoo{t] s) 

A simple computation shows that 

[C^{t),A{^t,Tpt)] = A{-/\^t + {V * \^t\^)vu ^Vt - {V * \^tf)^t) 
and therefore that 

-U:,{t;s)A{ipt,lpt)Uoo(.t;s) = 0. 

Hence 

^L{i;s)A{ipt,Tpt)Uoo{t;s) = A{ips,^s) 

and thus, from ()4.22p . 



□ 
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Next, we control the growth of the number of particles operator and of the kinetic energy operator 
with respect to the full fluctuation evolution U{t; s), defined in (j2.9p . 

Proposition 4.2. For every j G N, there exist constants C,K > (depending only on j, on D and 
on Wflln^ ) such that 

{U{t;s)i>,M^U{t-s)^) < Ce^^'-'^{i;,{Af+lf^+^ilj) (4.23) 

for all ip € J-, and all t, s G M. Moreover, there are constants C,K (depending only on D and on 
\\(p\\jj2 ) such that 

{U{t;s)i>,)CU{t;s)ij) < Ce^l*"'! (V', (/C + AA^ + 1) ^) (4.24) 

Proof. Eq. (j4.23p is taken from [2] [Proposition 5.1]. We show next the bound (j4.24p on the growth of 
the kinetic energy. To this end, we remark that, by Lemma [4.4l there exists a constant C (depending 
on D, W^pWh^) such that 



{U{t-s)il:XU{t-s)i)) <2{U{t-s)i),C{t)U{t-s)iP) + C{U{t;s)iP,{N+l) ^ + ]U{t- s)^) 



< 2 {Kit; s)^, C{t) U{t- s)iP) + Ce^l*-^l {iP, {M + 1) V> 



In the second inequality, we used (|4.23p . To bound the first term, we observe that 

^ {U{t; s)V, C{t)Uit; s)ij) = (u{t; s)V, C{t) U{t; s)^ 

where 



(4.25) 



C{t) = J dx {V* {ipt'ipt + (pm)) (x) ala^ 

+ j dxdyV{x-y) {iptix)Tp^{y) + ipt{x)Tpt{y)) a^Qy 

+ 2 dxdy V{x - y) {ipt(x)ipt{y) ala*y + Tp^{x)Tp^{y)axay) 



(4.26) 



+ / dxdy V{x - y)al {^t{y)a*y + (^t(y)ay) 



To bound the last term, we observe that, for any ip £ J^, 
^ dxdyV{x -y)(ft{y) {il;,ala*axi;) 

< ^= / dxdy \V{x - y)\ \(ptiy)\ \\ar,ay 



N 

< I dxdy\V{x - y)\\ipt{y)\'^\\axipf + — / dxdy\V{x-y)\\\ayacc^ 



< 



<c^^,(/c + caa(i + ^))v' 

for a constant C depending on ||(/?||j|^2. Here we used the bounds (|4.34p and, in the last step. 
The other terms on the r.h.s. of ()4.26p are just the time derivative Coo(t) of the generator Coo{t) of 
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the limiting dynamics. From the proof of Proposition I4.H we have Coo{t) < C{N + 1). Hence, we 
conclude that 




From Gronwall lemma, we find that 



\{U{t;s)^,C{t)U{t;s)ij)\ < Ce 



Inserting this estimate in (|4.25p . we obtain (|4.24p . 



□ 



Compared with Proposition 14.11 in Proposition 14.21 we do not have control on the growth of the 
square of the kinetic energy operator. Due to the cubic and quartic terms in the generator oiU{t; s), it 
is not clear whether a bound of the form (j4.6p holds with Uoo (t; s) replaced hyh({t;s). Fortunately, we 
only need some control on the growth of products of the form (AA+ iyiC{M+iy . For these operators, 
weak A^-dependent bounds are established in the next lemma going back to the original definition 
of the fluctuation evolution U{t\s) = W*{-/Nipt)e~'^'^^^^~^^W{-\fNips) (stronger, A^-independent, 
bounds for the growth of operators of the form [N + iyfC{M + 1)-^ with respect to the limiting 
dynamics Uooit^s) follow from Proposition 14.11 since, by Cauchy-Schwarz, {M + iyiC{J\f + 1)-' < 



{^M*{t; s){N + ly {JC + 1) (AA + iyu{t; s)ij) 
<C[iN + if^ {i',U*it; s) K,U{t- s)il^) + {iIj,{N + N + ly {K + + l) {N + N + ly V')] 



Proof. We proceed by induction over j G N. For j = 0, the claim is clear. Now we assume it to hold 
for j = n — 1, and we prove if for j = n. To this end, we observe that, from (j2.8p . 



where we recall the notation A{f,g) = a{f) + a*{g). Therefore, fixing for simplicity s = 0, we find 



Lemma 4.3. For every j G N, there exists a constant C such that 



W*{^ipt)MW[^ipt) = {M- ^A[ipt,Tp^) + N) 



(4.27) 
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By Cauchy-Schwarz, we have 

0) (M + ly (/C + 1) (AA + iyu{t; 0)V'^ 
< ^V, W* {VNip)e''^^^ MW{VN^t) {M + iy~^ (/C + 1) (A^ + l)-?"^ 

+ iV2 0) (A^ + ly-^ (/C + 1) (AA + 0)^) 

= 1 + 11 + 111 

In the last term, we can directly apply the induction assumption. We find 

{N + l)2{i-i) {'ilj,U{t; s)lCU{t- s)ij) 



III < CiV^ 



< c 



{N + lf^i^,U{t;s)}CU{t;s)tP) + {tP,{Af + N + iy (/C + AA^ + l) (AA + iV + l)J» 



.29) 



To bound the first term on the r.h.s. of (I4.28p . on the other hand, we use the fact that M commutes 
with TLn- From (jl2ZD> we find 



I = ^(AA + ^/NA{ip,Tp) + N)^P, U*{t- 0) {M + ly-^ {JC + 1) 

X (AA + ly-^ U{t; 0) [Af + VNA{ip, ^) + iV) 

< (AA^, 0) (AA + (/C + 1) (AA + 0) AA V> (4.30) 

+ iV (A(v9,^)^, U*{t; 0) (AA + l)^'-i (/C + 1) (AA + l)J'-i^/(t; 0) ^((/^, (?)V'> 
+ {ij, U*{t- 0) {N + 1)^'-^ (/C + 1) (AA + ly-^ U{t; 0) V') 

= Ai + A2 + A3 

The term A3 can be handled like the term III in (|4.29p . To bound Ai, we use the induction 
assumption. We find 

Ai < C 



(Ar + i)2(j-i) {Af^^U*{t;0))CU{t;0)Afi!) 

+ {Mi;, {M + N + iy~\lC+M^ + 1){M + N + ly-^Mtp) 
< C7(iV + l)2{J-i) (Mij,W*{VNip)e^'^^^ ICe-'^''^W{^/N^)J\fi;^ 
+ C{i;,{M + N + ly (/C + AA2 + 1)(AA + TV + l)^^') 

where we used that, by (12. 8|) . 

W*{VNipt)IC Wi^ipt) = (/C + VNA{Aipt, Alpt) + N\\Viptf) < /C + N\\Viptf <1C + CN. 



(4.31) 
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To bound the first term on the r.h.s. of (j4.3ip . we use Lemma |4.4[ We find 

< C {Mi), (/C + iV + 1 + N-'^M'^) Mi)) 

which implies that 

Ki<C {ij,{M + N + ly {K + M'^ + l) {M + N + ly i;) 

To estimate the term A2 on the r.h.s. of (j4.30p . we proceed similarly. From the induction assumption, 
we find 

A2 < C 



+ N {A{^, Tp)ij,{M + N + iy-\K: +M^ + i)iM + N + iy~^ a{^, Tp)i)) 

< C{N + (^A{^, ^)ij, W* {VN^)e^'^^*' {K + CN) e-''^^^W{^^) A{^, ^) i) 
+ CN {A{^,^)i), {M + N + iy-^{lC + M^ + 1){M + N + iy-^A{^,^)i)) 

Using Lemmas 14.51 and 14.41 we conclude that 

A2 < c {ip,{M + N + ly {K + M"^ + i){M + N + ly i)) 

Hence 

1<C\{N + {iPM{t; s)lCU{t- s)i)) + {i),{M + N + ly [f+M'^ + l) {M + N + ly i^) 
Finally, we estimate the term II on the r.h.s. of (|4.28p . We find 

II < N{i), VF*(\/iV(^) e^^'^* a*{^t) W{^^t) {M + l)-''"^ (/C + 1) (AA + ly-^ 

+ N(^i), H^*(^/iV(^)e*^^* a{^t) W{^ft) {M + iy~^ (/C + 1) (AA + ly-^ 

W*{^^t)a*{^t) e^^^* W{VNip)ij^ 

< N(^i)M*{t- 0)(a*((^0 + ViV) (AA + ly-^ (/C + 1) (AA + iy-\a{^t) + ^)U{t- 0)V^) 

+ N(^i:M*{t- 0) {a{^t) - Vn){M + ly-^ (/C + 1) (AA + iy-\a*i^t) - ^)U{t- 0)^^) 

< N^(^ij,U*{t; 0){M + ly-^ (/C + 1) (AA + ly-^ U{t] 0)V') 
+ N(^iP,U*{t; 0){M + 1)^'-^ (/C + 1) (AA + iyU{t; 0)^) 

where, in the last inequality, we used Lemma 14.51 By Cauchy-Schwarz, we find 

II < CN^(ij,u*{t; o){M + iy~^ (/c + 1) (AA + iy-^u{t; 0)^ 



(4.32) 



+ lU,U*it; 0)iM + 1)^ (/C + 1) (AA + iyU{t; 0)iP 
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The first term can be bounded as in (j4.29p . From (j4.28p we conclude that 

0) {M + ly (/c + 1) (AA + ly u{t- 0)?/^^ 
+ ]^I^^M*{t- 0) {N + ly {K + i){N + iyu{t- 0)^ 

This concludes the proof of the lemma. □ 

We conclude this section with two technical lemmas, which played an important role in the proof 
of Propositions 14.11 14.21 and of Lemma 14. 3[ The first result is used to compare the kinetic energy 
operator K, with the generators Hn-, and Coo{t) (these bounds are important, because it is much 
easier to bound the growth of the expectation of the generator than the growth of the kinetic energy). 

Lemma 4.4. Suppose that V'^{x) < D{1 — A) and let T-Ln be the Hamiltonian defined in \2. 5|) . C{t) 
the generator h2.12\l of the fluctuation dynamics U{t;s) defined in \2. 9^) and Coo{t) the generator 
\2.13\) of the limiting dynamics Uoo{t;s) defined in ^2.14^ . Then there exists a constant C > 
( depending only on D) such that 

K - C{M +l)<C^{t)<K + C{M + 1) , 
_/C -CAf[l + j^j<nN<2}C + CAT [l + ^ J and ^^33^ 

_/C - dM + 1) (^1 + —J < C{t) < 2/C + C{M + 1) (^1 + — 
Proof. The first bound in proven in [2] [Lemma 6.1]. To show the second estimate, let 

>V = ^y" dxdy\V{x - y)\a*a*yayax 



\Vix)\ <e{l-A) + 



From the bound V'^{x) < D(l — A), we obtain that, for every e > 0, there exists > with 

1 

e 

Therefore, the restriction W*-"-* of W on the n-particle sector is bounded by 

Kj j = l 



Choosing e = 5N/n, we find 
Hence, as operators on the Fock-space, 



W<5{}C+Af) + ^ (4.34) 



With 5 = 1/2, we conclude that 



1 1 7V^ 
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The upper bound for can be proven analogously. To prove the last estimate in (|4.33p . we observe 
that, by Cauchy-Schwarz, 

C{t) = Coo{t) + j dxdyV{x - y) a* {a*y(ft{y) + ayTpt{y)) "-x + j dxdy V{x - y) ala^aya^ 



1 



>Coo{t)-2 j dxdy\V {x - y)\\Lpt{y)\ ala^ - — J dxdy - y)| a* 0*0^03; 

The lower bound for £(t) follows now by the lower bound for Coo{t), by the bound ()4.34p . and since 
clearly 

J dxdy\V{x-y)\\^t{y)\^ala,<\\\V\*\^t\^\UAr <CAf 

The upper bound for C{t) follows similarly. □ 

Finally, in the next lemma we show how to bound creation and annihilation operators by the 
number of particles operator after commuting them through the kinetic energy. 

Lemma 4.5. For j G N, there exists a constant C , such that 

a*{^){M + iy]C{M + iya{ip) < ||v?fAA-''+i/2/cAAJ'+i/2 

a{ip){M + iyjc{M + iya*{ip) < c \\ip\\l^ {M + iy+^/\ic + i){M + ly+^z^ (^-^^^ 

Proof. We observe that, if /C^") denotes the restriction of IC onto the n particle sector, 
i;, a*i^){M + lyiCiM + ly a((^)0) = ^ ^ ((AA + l)^a((^)V') , /C^^^ {{M + iya{ip)ij) 

n>0 



n „ 

^{n + lf^Yl J dxi...dxn |v^^.(a(^)V')(")(xi,...,x„) 



n>0 j=l 

Since 

(a(v.)7^)(") {xi,...,Xn) = {n + 1)1/2 I ^(n+l) (a,^ a,^^ . . . ^ 3,^) 

we conclude that 

n „ 

= ^^{n + lf^^^^^ j dxi...dxndxdyi^{x)ip{y)V^^ilj^"-^'^\x,xi,...,Xn)VxjilJ^'^'^^\y,xi,...,Xn) 

n>0 j=l 

n „ 

< Y.^n + lf^+^Y. / dx^---dXndxdy\^(x)\''\V,^i^^^+^\y,x^,...,Xn)\'' 

n>0 i=l 

n+1 „ 

< M?^{n+lf^^^Y^ I dxi...dXndXn+l\yx,i^^'^^^\xi,...,Xn,Xn+l)\'^ 

n>0 i=l 

To prove the second bound in (j4.35p . we commute a*[(p) to the left and 0(93) to the right of fC. Since 

a{ip)K.a* {if) = a{ip)a* {ip)K, + a{ip)a*{—/S.ip) 

= a* {(p)a{(f)IC + \\(p\\'^ IC + a{(p)a* {— A(p) 

= a* {(p)}Ca{ip) +IC + a*{ip)a{-Aip) + a{ip)a* {-Aip) 

= a*{ip)lCa{^) + yW^JC + a*((^)a(-Av9) + a*(-Av9)a(v9) + 
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we find 



= {M + 2ya{ip)}Ca* {ip){M + 2y 

<{M + 2ya*{^)ICai^)iM + 2y + \\ipf {M + 2yiciM + 2y 

+ {M + 2y ia*{^)ai-A^) + a*{-A^)a{ip)) {M + 2y + \\Vipf{M + 2f^ 

< a*{^){N + SyiCiM + Syaiip) + (AA + 2y {yf K, + ||V(^f ) {N + 2y 
+ {M + 2y (a*((^)a(-A(^) + a* {- A^)a{ip)) {N + 2y 

Using the first bound in (j4.35p to control the first term, and the observation that 

{N + 2y {a*{^)a{-A^) + a*(-A(^)a(vp)) [N + 2)>) 

< 2||a((^)(AA + 2y^\\ \\a{-A^){M + 2^ V^|| 

< Wvfi^, {N + 2)2^'+V) + llVv^f {M + 2yK{M + 2)H) 

where in the last inequality we integrated by parts to move one of the derivatives in — A(/p back onto 
ip. We conclude that 

a{^){M + iyic{M + iya*{ip) < c\\ip\\l^ {Af + iy+^/^{}c + i){Af + iy+^/^ 

□ 



5 Regularization of the potential 

As in [2], our analysis requires a regularization of the interaction potential. This is needed to establish 
the convergence of the fluctuation dynamics U{t;s) to its limit Uoo{t]s). We introduce therefore an 
A'^-dependent cutoff in V , vanishing sufficiently fast in the limit N ^ oo. The last condition should 
guarantee that the quantities we are interested in do not change when V is replaced by the regularized 
potential. On the other hand, the presence of the small cutoff will give us an (A^-dependent) L°°- 
bound on the potential which is crucial in the analysis (in particular, in Section [6|). 
For an arbitrary sequence > 0, we set 

V{x) = sgn(F(x)) • min{|y(x)|,a^^} (5.1) 

where sgn{V{x)) denotes the sign of V{x). Note that, by definition |V^(a;)| < . Moreover, the 
assumption V'^{x) < D{1 — A) implies that the same bound holds also for V, independently of A^: 

V\x) < D{1 - A,) . (5.2) 

We define the regularized Hamiltonian 

N ^ N 

j=l i<j 

On the Fock space J^, we define 

Hn = J dxVxalVxax + ^ ^ ^^^2/ ^(^ ~ v) ala^Uyax (5.4) 
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and the fluctuation dynamics 

U{t- s) = e-^'^(*'") W*{^/N^t)e-'^''^*-'^W{^/N^s) (5.5) 

with the phase 

^{i'^s) = — J dr J dx{V *\(pr\'^){x)\(pr{x)\'^ 
and where ipt solves the regularized (A^ dependent) Hartree equation 

idtift = -Aipt + (V * l^tH ipt (5.6) 



We denote the generator of U{t;s) by C{t). Clearly, C{t) has the form (12.121) . but with V and (pt 
replaced by V and, respectively, ^pt- We also modify the observable Ot defined in (|3.ip by setting 



N , N 



^^■■=^Y. [O^'^ - = E (O^^'^ - On)) (5.7) 



To compare the original many body evolution (generated by the Hamiltonian Hi\}) with the 
regularized many-body evolution (generated by -f/^iv), we use the following lemma. 

Lemma 5.1. Let ipN,t = e'^-^^V®^ a^d ijN,t = e-*^^V®^> with 99 G H^{m^) such that = 1. 
Then there exists a universal constant C > such that 



<CNaN\t\ (5.5 



for allN e'N,te 



The proof of this lemma can be found in [2] [Lemma 2.1]. We can also compare the solution of 
the original Hartree equation ()1.3p . with the solution of the regularized Hartree equation (15. 6p . 



Lemma 5.2. Let if £ H^{W^), with Hvlb = 1- Let ift and ^pt denote the solutions of the Hartree 
equations U.3\) and, respectively, \5. 0(1 . with initial data (/?. Then there are constants C,K > such 
that 

\\n -^t\\< CaAre-^l*! 

for every t G M. 

The proof of this lemma can be found in [2] [Lemma 2.2]. 

As a consequence of the last two lemmas, it is enough to show (j3.4p for the regularized quantities, 
assuming that the cutoff is sufficiently small. In other words. Lemma [3.21 follows by the next lemma, 
where the original potential V is replaced by its regularized version V. 

Lemma 5.3. Fix k £N. Choose an integer r > 3 + k/2 and let un = . Let V , Hn, 'ft O'^d Ot 
he defined as in (EJ]), (ED, {SJ), fJ^j respectively. Let i)N,t = e'^^'^^ip®^ . Then 

~, / ~ ~ ~ \ I if k is odd 

lim IE,^„ Of = lim UN,t,O^^P^A = kl (5.9) 



AT^oo >iv.-* - ^'^'^ V^^'^^ ^i^.t, - I J^af if k IS even 

with at as defined in 113.^) . 

Using Lemma 15.31 we can prove Lemma 13.21 (and thus Theorem 13. ip as follows. 
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Proof of LemmalEM Fix G N, t G M. Choose r > k/2 + 3 and let oat = iV Observe that 



\\Ot\ 



1 ^ 



< ^ ||0 - E<p,0|| < 2Vn\\0\ 



The same bound holds for Ot- Moreover, from Lemma 15.21 we have 

\\0t-dt\\<^\{ipt,0^t)-{^t,0^t)\<2^\\ipt-M\\0\\ <CVNaNe^\'\ 
This implies that 

for an appropriate /c-dependent constant C. Therefore, using Lemma |5.H we find 



< 



< 



+ 2 



'4'N,t - '>pN,t 



\Ot\ 



which converges to zero, as N ^ oo. We conclude that 



and therefore 



N~>-oo 

follows from Lemma 



lim nf = lim Er Of 



□ 



6 Comparison of Dynamics 

Finally, we will need to compare the (regularized) fluctuation dynamics U{t; s) with the limiting 
dynamics Uoo{t',s) defined in ()2.14p (note that the limiting dynamics is defined in terms of the 
original potential V, and the solution of the Hartree equation (jl.3p . without cutoff; the reason is 
that, in the limit — )• oo, the A^-dependent cutoff oat disappears). 

Proposition 6.1. Suppose that, in the definition (5. 1\) of the regularized potential V, the cutoff a n 
is such that > A^"'', for some r G N, r > 1. Then, for any j G N/2, there exist constants 
C,K > depending on r and j, such that 



(M+iy [U{t;s)-U^it;s))^ 

< C7e^l*-^l + aN{N + iy+^ 



^^_^{4r-2)(4+2i)+6^1 ' (AA + 1 



1/2 



\2j+2 



(6.1) 



Note that the bound (16. ip is only useful for j < r — 1; it will be applied under this condition. 

Proof. We fix t > and s = (all other cases can be treated analogously). We use 

ft 



U{t;^)-U^{t;Q)= 1^ dTUoo{t;T){/:{T)-Coo{r))U{T;0) 



(6.2) 



26 



where C{t) is the generator of the regularized fluctuation dynamics U{t]s) and has the form (j2.12p . 
with V and ipt replaced by V and ipt- Hence 



(AA+l)MZ^(t;0)-Z^oo(t;0)U 



t 

< / dr 

•t 







dr 







where we decomposed 
with 



dT 



C{t) - £oo(r) = £2(r) + £3^ + £4 



C2{r) = Jdx(v* \^t\\x) - V * \ipt\\x)^ ala, 

+ j dxdy {y{x - y)ipt{x)if{y) - V{x - y)(pt{x)Tp^{y)^ a*ay 

+ / dxdy (V{x - y)^t{x)^t{y) - V{x - y)ipt{x)ipt{y) \ a*a* + h.c 



and 



>C3(t) 



1 



dxdyV{x - y)al [ft{y)a*y + ft{y)ay\ a 



£4^ = — J dxdy V{x - y) ala^ayax 
We begin by estimating the first term on the r.h.s. of (j6.3p . From Proposition 14. II we have 



{N +iyU^{t-T)C2{T)U{T-Q)ij 



[N +iyC2{T)U{T-Q)ij 



We write 
where 



^1, 



A 



2,T 



dx{y*\^t?{x)-V*\^t?{x)) ala, 

dxdy {V{x - y)^t{x)^{y) - V{x - y)ipt{x)Tp^{y)] a^a. 



with 



^3,r = Br + B* 

dxdy {v{x - y)^t{x)^t{y) - V{x - y)ipt{x)^pt{y)^ ala*y 



B. 



From ()6.4|) . we find 

{N + iyUoo{t;T)C2{T)U{T-Q)ij 



Ai^r iM + iyU{T;0)7p 



Br {M + 2,yU{T- O)'0 



A2,r {M + iyU{T- 0)^- 

Bl {M - iyU{T- 0)i; 



(6.3) 



(6.4) 



(6.5) 
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Using \V -V\< aN\V\'^, the assumption < ^^(l - A), the propagation of regularity bound ()4.3 
and the bound \\Lpt — ^t\\2 < CaAre^'*' from Lemma 15.21 we find a constant C such that 



sup 

X 



{V*\^r\^){x)-{V*Wr\^){x) 



This bounds imphes that 

Similarly, we find A^^ < Ca'j^e^'^ Af"^ . To bound the last two terms on the r.h.s. of (j6.5p . we write 



Bt= J dxdy (V{x -y) - V{x - y)^ ipt{x)ipt{y) <a* 

+ J dxdy V{x - y) (^t(x) - <ptix)) (pt{y) ala*y + j dxdy V{x - y)'ft{x) {^t{y) - ^t{y)) ala^ 
Cauchy-Schwarz inequality implies 

(■0, BrB^i^) < 1 dxdydzdw (v{x - y) - V{x - y)) (y{z - w) - V{z - w)) 



+ j dxdydzdwV{x - y)V{z - w){ipt{x) - ipt{x)){(p^{z) -Tpf-{z)) 
+ / dxdydzdwV{x - y)V{z - w){^t{y) - Vt{y)){'?t{w) -Wt{w)) 



(6.6) 



X Lpt{x)ipt{z) {il!,ala*a;,ayjip) 
= 1 + 11 + 111 

Using |y — 1^1 < unV^ and the operator inequality V'^ < D{1 — A), the first term is bounded by 
l<a% dxdydzdwV'^{x - y)V'^{z - w)\(pt{x)\\(pt{y)\\(pt{z)\\(pt{w)\\\ 
<a%j dxdydzdwV"^ {x - y) ^^(z - (z)p|^t('«^)P 
<Ca% j dxdy xa-xayipW"^ + WaxayipW"^) 

The second term on the r.h.s. of (j6.6p . on the other hand, can be estimated by 



II < y dxdydzdw V"^ {x - y)\Lpt{y)\'^\'^t{x) - ipt{x)\'^ \\ 

where we used Lemma |5.2[ The third term on the r.h.s. of ()6.6p can be handled similarly. Hence 
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Analogously, one can prove that 

(the factor 1 takes into account the contribution of the commutator between Br and B*). We 
conclude from (16.51) that 



{N +iyUoo{t;T)C2{T)U{T-Q)ij 



< CaN e 



Kt 



+ Cajy e 



Kt 



<CaN e^*|| (AA+1)2J+3V;|| 



(6.7) 



Kt 



< CuN {N + ly+^e^' (/C + + 1)1/2 ^ l)2i+2^ 

where we applied Proposition 14.21 and Lemma 14.31 (and where \j + 1/2] denotes the smallest integer 
larger or equal to j + 1/2). 

To bound the second term on the r.h.s. of (j6.3p . we use the inequality 

£3(t) (AA + if^C^it) < ^(AA + 1)2^+3 
which is proven in [2j [Lemma 6.3]. This implies that 



(AA+l)^^/oo(t;r)£3(r)^/(r;0)V 



< C7Ar-i/2gX(t-r) ^ iy+^/^U{T; ov- 

Finally, to bound the last term on the r.h.s. of (|6.3p . we note that, assuming > 
{Af + lycliAf + ly < CN-^ f/C + AA2"(2i+4) ^ 



(6.8) 



(6.9) 



This is in fact the only place in the proof where a cutoff in the potential is needed. To prove (j6.9p 
we notice that, when restricted to the n-particle sector 

1 



(n)- 



Ar2 



< 



[n + l)2j+3 1 (n < iV2j+3 



Ar2 



(n + l)2j+4i U > 



aim 



a 



N 



Assuming > N we find 

(AA + l)2i+4i(AA> A^W) 



A2a2 



< A 



2r-2- 



27+3 (AA + l)2i+4+^ 



29 



for any £ > 0. Choosing I = (2r — l){2j + 3), we get (|6.9p . This imphes that the third term on the 
r.h.s. of (|6.3p is bounded by 



1/2 



^^^4r(4+2i)+2^-^ ^ 



1/2 



(6.10) 



where we used again Proposition 14.21 Inserting (j6.7p , (|6.8p , (|6.10p in (j6.3p , we obtain (|6.ip . □ 



7 Computation of the moments 

In this section we prove Lemma [5.31 To this end, we fix /c S N and we choose an integer r > A;/2 + 3. 
We let Oat = A^"*", and we define the regularized potential V as in (jS.ip in terms of oat. For the 
given compact operator O on L^(M^) we introduce the notation Oc = O — {tpt,0^pt)- Later, we will 
also use the notation Oc = O — {ipt, Oipt) (as a rule, notations with a tilde denote quantities defined 
in terms of the regularized potential V and the solution ipt of the corresponding Hartree equation 
(lESD). We write 



V'JV.t ' 



1 

Tk/2 ^ 



N 



]S[k/2 ^ , ^ 

m=l ^ ^ iiH |-«m=fc 



E 



^l' • • • ^m' 



The sum over ii, . . . ,im runs over all integers ii, . . . , im > 1 with zi + • • • + = /c. We separate next 
now the indices ii, . . . ,im which are equal to one (r denotes the number of such indices). We write 



WN,t * 



k\ 



where the sum runs over ji, . . . , jm-r runs over all integers ji , . . . , jm-r > 2 with • ■+jm-r = fc— 
In the following, we denote by ^pN,t the vector {0, . . . , 0, 11)^,1-,^^ • • • } in the Fock space F. Then, if 
O* (x; y) denotes the integral kernel of the operator O* , we obtain 



E~ Ol 

WN,t * 



1 ^ 
Tk/2 ^ 



1 / TV \ / m 



E 



jVfc/2 ^ /V™ V ^ V y ^ Jl! • • • jm-r! 



a*(^)^ 



rj, e^^'^* a:, . . . a* o^;' ... e~^^^^ ^^ll n 



with the (regularized) Hamiltonian Hn defined in ()5.4p . Next, we notice that 

a%)f 



-Jl = dNPNW{VN^)n 
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where W{VlV(p) is the Weyl operator defined in (j2.7p . Pn is the orthogonal projection onto the 
iV-particle sector of the Fock space and dN = \/Me^/2j\r-Af/2 ^ ^i/4_ Denoting by U{t;0) the 
(regularized) fluctuation dynamics introduced in ()5.5p . we find 



m=l / r=l ^ ^ 



J l ! • • • jm—r ! 



X y ds^idx ]^ ... dXm'j^S^m Oc(3^li ^l) ••• Oc(2^r) 2Jr)^c^ -^r+l) ••• ^c^^ '^(■^mj^^m) 

X 0) «^ + VN^ti^i)) . . . «,„ + VN^tixm)) 



where we introduced the notation 



Note that, from [2^^, Lemma 7.1], we have the estimate 

||(AA + l)-i/2^^|| <C7 



(7.11) 
(7.12) 

(7.13) 



for a constant C > 0, uniformly in N. When we expand the product of the 2m parenthesis on the 
r.h.s. of (17. lip , we find several terms, with different numbers of factors ipt or if^. Each term has the 
form (after appropriate change of variables, and with a constant C depending on k, m, r,ji,... ,jm~r) 



CN 



a + b+2d-k 



N' 



m 



{Lpt,Ol^ipt) ■ ■ ■ {ipt,Ol'^ipt) / dxidx'i . . .dxndx'nOl^{xi,x'i) . . . 0^"(a;„,a;^) 



X {U{t- 0)Civ, a*(0,'Vt) • • • a*{Ol-^t)al, . . . a^^a,,^ . . . a,, a(0,^Vt) • • • a{Ol^^t)U{t- 0)17 



(7.14) 



where n,a,b,d are integers with n + a + b + d = m, and {pe}i^i, {Q£}e=i,{f'£}i=i, {s^}^=i are 
appropriate sequences of integers (r of them are equal to 1) . Note that a + 6 + 2d is the total number 
of factors cpt and ; since each such term carries a weight N^^"^ , this explains the appearance of the 
factor Ar(a+fe+2<i)/2)_ r^^^ absolute value of (fTlip is bounded above by 

CN'^'^Wdcir I dx[ . . . dx'^\(l({t;0)CN,a*{d^,'ipt) . . . a*{dr^t) 

X a*(0?^ (., x[)) . . . a*(0?"(., x;))a,. . . . a,, a(O^V*) • • • a{dl^^t)m 0)0 
< CN"^"^ ||Ocir''^+-+^^'* j dx[ ...dx'^ 

{M + l)2+(n+-)/2 . . . aidl^ipt) . . . a{dl^^t)U{t- 0)Q 
(AA+ l)-2-{«+-)/2o(OPn(.,4)) . . . a{dl-{.,x'^))a(dl^^t) . . . a{d^,^^tMt-0)^N 

Q^||2(si + ...+s,)+2(n+...+rd) II (AA + l)2+"+('*+^)/2z7(t; 0)1711 



a + b+2d-k 

< CN 2 



+ CiV^^^^||0,|pto+-+«") / dx[...dx'JOP"i.,x'jf ...\\OP^{.,x[)fm + l)-'U{fMNf 



a + b+2d-k 

< CN 2 



1^ ||2(siH hst+riH hr-rf) , ii^ ||2(cjiH l-9a) ||/-)Pn ||2 

I'-^cll -r IIL/cll W^c llHS 
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where in the last hne, we used Proposition l4.2l (which imphes that ||(AA+1) '^U{t] 0)(A/'+l)^/^|| < C). 
From (j7.13p . the absolute value of (j7.14p is thus controlled by 



„ a + b+2d-k 

CN 2 



^||Q||2(siH hS(,+riH f-rd) _|_ ||Q||2((jriH h^a+PiH Vpn-n) ||(9||^^ 



Hence, we conclude that only the terms with k < a + h + 2d give an important contribution in the 
limit N — )• oo. This implies that m > k/2 (and m > (A; + l)/2, if k is odd). Moreover, terms where 
both <ft{xk) and ^pt{xk) appear, for some k = l,...,r, vanish because {(pt.,Oc^Pt) = (this is the 
reason for separating the indices which are equal to one). These two observations imply that only 
terms with r < 2m — k may contribute in the limit — )• oo (if r > 2m — A; + 1, terms with at least 
k factors of and ^p^ will contain at least one factor {ipt^Oc'f't) and therefore will vanish). On the 
other hand, since j'l, . . . ,jm-r > 2, we have 

k-r = h jm-r > 2{m - r) 

which implies that r > 2m — k. We conclude that only terms with r = 2m — k, with Ji = • • • = 
jm-r = 2, and with exactly k factors of (ft can contribute in the limit N ^ oo. More precisely, after 
integration, we obtain 



V'JV.t ' 



k 

y — 



m=k/2 



+ o(l) 



m 



m\\2m-k 2^=-"^ 



k\ 



k—m 



X (e7v,Z^*(t;0) : {a{0,^t) + a*{0,^t, 



2m— k 



■.U{t;0)n) 



(7.15) 



as iV — )• oo. Here : • : denotes the operation of normal ordering of the creation and annihila- 
tion operators; all creation operators have to be written on the left, and all annihilation operators 
have to be written on the right, as if they commuted (for example, : a(/i)a*(/2)a*(/3)a(/4) : = 
'2*(/2)a*(/3)a(/i)a(/4))- Using Proposition 16.11 and (|7.13p . we can replace, in (|7.15p . the (regular- 
ized) fluctuation dynamics U{t;0) with its limit Uoo{t',0), defined in (I2.14p . In fact. 



^N,l^*{t;0) : [a{Oc^t) + a*{Oc^t)) 



2m— k 



2m— k 



■.l(oo{t;0)n 



< 



2m— k 



■.u^{t;o)n 



+ 



2m— k 



U{t;0) -Uoo{t;0))n 



2m+k 

<\m+i)-'/'(N\\ 

(AA + 1)1/2 ^i{*(^t;0) -UUt;0)) a*{d,^tya{d,iptf^-''-^U^{t;0)n 
(AA+ l)i/2^(t;0) a* (dc^tYaidc^tf {u{t;Q) -U^{t-Q)) Q 



(7.16) 



2m — k 
3 
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To bound the first term in the parenthesis, we use Proposition I6.1i We find 

< Ce^*(iV-i/2 + a^(Ar + i)2) 

X I (/C + + 1)1/2 ^ lfa*{dc^^ya{dc^t?"'-''-moo{t; m 

Using Lemma 14.51 we conclude that 

< Ce^*(iV-i/2 + ajv(iV + 1)^) + -^^^'"^ + l)^^^!-^ + l)^+'""''/^^oo(t; 0)17 

|/CZ^oo(t;0)17|| + (AA + l)6+2™-fc^/^(t;0)J^ + (AA + l)iO'^+i+'"-'=/2^/^(t;0)n 



where, in the last inequality we used Cauchy-Schwarz to separate K^^"^ and {J\f + 1)^+™ ^Z"^. From 
Proposition 14. 11 we obtain 



1)1/2 fK*^t;0) -W^it;0))a*{Oc^tya{Oc^t) 



\2m--k- 



~ \\2m~k 



< Ce^*(iV-i/2 + ajv(iV + 1)2) \\0, 

< Ce^*(iV-i/2 + ajv(iV + 1)2) (||V 0(1 - A)-i/2|| + \\0 



2m— fe 



Since on = N for r > 3 + k/2, the r.h.s. of the last equation converges to zero, as — )• oo. To 
estimate the second term in (j7.16p . we observe that, by Proposition 14.21 and Proposition 16.11 

(AA+ l)i/2ZY(t;0)a*(Oe^t)^a(Oe^t)2— (u{t;0) -Uoo{t;0)) O 



< Ce 



{M + If a*{Oc^tya{Oc^tf"'-''-^ {U{t-0) -l(oo{t;0)) n 
||(AA+1)2+— '=/2 (U{t;0)-U^{t;0)^n 
< Ce^* (iV-i/2 + aN{N + lf+"^-k/2^ \\o 



< Ce^*||Oc^t||2"^-'= 



1 2m— k 



which converges to zero as — )• oo, because = N ^, with r > 3 + k/2 > 3 + m — k/2. From 
(j7.16p . using also (j7.13p . we conclude that 



IpN* ' 



k 

y — 



m=k/2 



+ 0(1) 



m 



! \2m - k 2^ 



k\ 



k—m 



2m— k 



(7.17) 



as — )• oo. 

Suppose now that /c = 2£+ 1 is odd. Expanding : {a{Oc'Pt) +«*(Oc(/7t))2™'~2^-i ; -^g g^d a sum of 
(normally ordered) terms of the form {a*{Oc'^t)Y{o-{Oc'-Pt)Y with p + q = 2m — 2i — 1. Conjugating 
with the limiting dynamics Uoo{t; 0), and using that 

U*^{t; 0)a{O,ipt)Uoo{t; 0) = a{Ft) + a*{Gt) 
U:,{t;0)a*{O,^t)Uoo{t;0) = a*{Ft) + a{Gt) 
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with Ft = U{t; 0)Ocipt and Gt = JV{t; O)0cipt, we find that K^{t; 0) : (a(OcV?t) + a* {Oc^t)f"''^^~^ : 
Uoo{t\ 0) is a hnear combination of terms of the form {a* {Ft) + a{Gt)Y{a{Ft) + a*{Gt)Y with p + q = 
2m — 2£— 1. Expanding the parenthesis, we conclude that it is a hnear combination of terms given by 
(not necessarily normally ordered) products of creation operators a* (Ft) and a*{Gt) and annihilation 
operators a{Ft) and a{Gt). In every term the total number of annihilation and creation operators is 
2m — 2i — 1. Consider a single term, with a total of p creation operators, and of ^ = 2m — 2i — l—p 
annihilation operators. Since these operators act on the vacuum, this term is a vector in the sector 
of F with exactly p — q = 2p — 2m + 2^+1 particles. On the other hand, from Proposition 18.11 we 
observe that the odd components of the vector = . . .} £ F vanish in the limit — ?■ oo. 

Hence, we conclude that, for any fixed £ G N, 

E- of +^ ^ 

■>PN,t * 

as A?" — )• OO. 

Next, we consider the even moments. For k = 2£, (|7.17p gives (changing the summation variable 
to m — £) 

f 

~ 2£} I 

^^.,^''= E 2m\l - m\2^-^ QcVt)''" (^^oo(t; : iaiO^^t) + a{p,^t)f^ : U^it- 0)0 

m=0 



+ o(l) 



as — )• oo. From Proposition 18.31 we find 

^^^pf = (Uoo{t; 0)eiv, {a{Oc^t) + a{Ocipt)f^Uooit; 0)1^) + o(l) 
With the notation A{f,g) = a{f) + a*{'g), and using Theorem 12. 2^ we find 
E^^_^Of = (u^{t;0)CN,A{Oc^t,JOcipt?'Uoo{t;0)n)+o{l) 

21, 



(7.18) 



^^TV, A {Uit; O)0cipt + JV{t- Q)OcH^t, V{t; Q)0,^t + JU{t; 0)Oc^tY' ^ 
From Proposition 18.11 we write 

j>0 



(7.19) 



Using LemmaE2]to expand A {U{t; O)0cipt + JV{t; O)0cipt, V{t; 0)Ocipt + JU{t; 0)Oc^Ptf^' ^, we find 

2£ £ 
j=0 m=0 

X {a*{ipyn, a*{U{t; O)0,(^t + JV{t- {))OcVt?"'^) 
Since, for any f,gG L^(M^), and for any j, G N, 

{a*{fyn,a*{g)''n)=6,k{f,g)''k\ 
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we conclude that 



m=0 



2^-™(f -m)! 

m=0 ^ ' 

From Proposition 18.11 we find that, for any m = 0, . . . ,£ fixed, 

A2m) ^ (-ir 

^iv 2™m! 

as — )• oo. Hence 
hm Er O?^ 

= m E(-ir(0 llf^(*;o)a¥'t + -/^(t;o)Oc¥'tf^'~'"^</',t/(t;o)Oc(^t + jm</'t)'" (7.20) 

= ^ {\\U{t; 0)Ocn + JV{t; O)0c^tf - U{t; O)0c^t + JV{t; O)0c^t)^f^ 
Now, we observe that 

-i {{if, -^), Q{t- 0) (0,(^t, JOc^t)) = - U{t; 0)Ocipt + JV{t; 0)Ocipt) 

+ V{t; 0)Ocipt + JU{t; 0)Ocn) (7.21) 
=2Im U{t; O)0cn + JV{t; 0)Ocn) 

From (f2J9]) . and since e(t;0)* = SBt'^S (with 5 as defined in dTTUIl l. we conclude that 
2Im ((^, ?7(t; O)0,(^t + jy(t; 0)OeV?t) = - i (e(t; 0)*(^, -^), {Oc^t, JOc^t)) 

= -i(e(t;0)-i((/p,(^),(O,v9t,-JOeV't)> (7.22) 
= - i {{ipt,Tp^), {Oc^pt, -JOc^t)) = 2Im {ift,Oc^t) = 

Thus 

and, from (I7:20]) . 

^lirn^E^^^^ Of = ^ (||^(t;0)Oe9^i + JV{t-^)0,iptf - \{^,U{t-^)0,^t + JV {t;Q)0,^t)?f' 

(7.23) 



with 

(j^t = \\U{t-0)Ocipt + JV{t-0)O,^t\\^ - \{^,U{t;0)O,vt + JV{t;0)OcVt)\^ 
Finally, we note that 

{{ip,lf),e{t;0){Oc^t,JOc^t)) 

= {i^,lp), {U{t; 0)Oc^t + JV{t; 0)OcVt, JU{t; O)0c^t + V{t; 0)Ocipt)) 
= 2Re {if, U{t; O)0cipt + JV{t; O)0cipt) 
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and that, similarly, 



\\U{t; 0)Oc^t + JV{t; 0)Oe(^tf = ^ 0) {OcVt, JOcVt) , @{t; 0) {Oc^t, JOc^t)) 



This gives 



<i = \ 0) {Oc^u JOcipt) , @it; 0) {Oc^t, Ocipt)) 

e(t;0) {OcVt,JOcn),^{ip,^)'^ (^-^{^,ip),e{t;0){Oc^t,JOcn: 



(7.24) 



Writing Oc = O — IE(^jO = O — {ft, 0(pt), we easily find 

^It = \ [(©(*; 0) {On, JO ft) , e{t; 0) {o^t, JO^t)) 

e{t; 0) {Oft, JOft) , ^(v''^)) o) {Oft, JOft] 

\\U{t-0)Oipt + JV{t-0)Oiptf - \{ip,U{t;0)Oipt + JF(t; 0)O(/Jt)|2 
a? 



(7.25) 



With (|7.23p , this concludes the proof of Lemma 15.31 



8 Combinatorial estimates 



In this section, we collect some technical results, used in Section [21 to handle the combinatorial 
problems. 

The next proposition is used to determine the properties of the Fock vector defined in (|7.12p 
in the limit of large A^. We use here an important observation from 



Proposition 8.1. Let 



with dN = e^/2^]vT Then 



with 



-£/2 



m 



j=0 



N-e + jie-jij\ 



Moreover, for any fixed m G N, we have ^l^'""''^'' — )• and 



(8.2) 



J2m) ^ (-1) 



2'"m! 



as N ^ oo. 
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Proof. We start from the definition of ^n, given by 



with djq = e^/2\/MiV"^/2 ~ iVi/4, Using the properties of the Weyl operator W*{^ip), we obtain 

^^N\ ^/M ^ j! 

The restriction of on the ^-particle sector is given by 



-£/2_ 



iV! 



:^) 



N-t + j\i-j\j\ 



a*{ipYn 



This proves (jS.ip . Now we study the asymptotics of as — )• oo. To this end we recognize, 
following [23], that 

eW=iV-^/2Lf-^)(iV) 



where 



j=0 



J J J! 



are the generalized Laguerre polynomials of dergree i, which satisfy the recursion relations (see, for 
example, [23]) 



a + l- X ^(a+l) 



This gives the recursion 



Ai) _ 1^ ^-1/2 .(^-1) _ 1 .(£-2) 



(8.3) 



A simple computation shows that = 1 and ^^^^ = 0. Eq. (j8.3p implies therefore that, for any 
fixed £ G N, Cff ^ if ^ is odd, and 



if £ is even. 

The next lemma is used to expand 



(-1)' 



(-1)' 



£(^-2)... 2 2^/2(^/2)! 



□ 



A {U{t; O)0c^t + JV{t; 0)0^9^*, V{t; 0)0^9^* + JU{t; O)0,iptf O 



{a*{U{t; 0)Ocipt + JV{t; O)0,(^t) + a{U{t; 0)O,(pt + JV{t; O)0c^t)f^ n 



appearing in (j7.19p . 
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Lemma 8.2. For any F S L^(R^), we have 

2m\e - m\2^ 



2/1 

(a(F) + a*{F)f^n = V , , \\Ff^^-"'^a*iFf"'n 

m=0 



Proof. We proceed by induction over /. The statement is clearly correct for / = 0. Assume it is 
correct for a fixed / G N. We prove it holds true for i+ 1. To this end we write, using the induction 
assumption, 

(o(F) + a*{F)f^+^n = (a(F) + a*{F)f (a(F) + a*{F)f^n 

= y \\Ff('-^''\a{F) + a*{F)fa*{Ff^n ^^'^^ 

m=0 

Next, we observe that, using the canonical commutation relations, 
{a{F)+a*{F)fa*{Ff'^ = a* {Ff'^"'+^^ + a{Ffa* {Ff"^ + 2a{F)a* {Ff'''+^ - \\Ffa*{Ff'^ 

= a*{Ff^"'+^'^ + 2m\\Ff a{F)a* {Ff""-^ + (4m + l)\\Ffa*{Ff'^ 
= a*(F)2('"+i) + 2m(2m - l)\\Ff a*(F)2('"-i) + (4m + a*(F)2'" 

Inserting in (j8.4p . we find 

(a(F) +a*(F))2^+2^ = V (^-"^)a*(F)^("+^)l^ 
V V ; V 2m!/ - m!2^-'" " " ^ ^ 

m=0 

^ 2/!(4m + 1) II p||2(^+l-m) */7^N2mn 

+ Z^2m!/-m!2^-™" " 

m=0 

+ Z^2(m-l)!/-m!2^-'»" " 

m=l 

Changing the summation index to m + 1 in the first and, respectively, to m — 1 in the third term on 
the r.h.s. of the last equation, we obtain 



{a{F) + a*{F)f^+^n 

^2(m-l)!/ + l-m!2^+i-™" " ^ ' 



m=l 



m=0 

^-1 



2m!/ - m!2^ 
2/! 

2m!/ - 1 - m!2^ 



+ V — llFlP(^+i-™)f7*('Fl2"0 



m=0 

o*(F)2(^+i)f] + (/(2/ - 1) + 4/ + l)\\Ffa*{Ff^n 

I V ||^||2(^+l-m) 

^ 2m!/ + 1 - m!2^+i-™ " " 

m=l 

X (2m(2m - 1) + 2(4m + 1)(/ + 1 - m) + (/ + 1 - m)(/ - m)) a*(F)2™0 



2(^ + 1)! ||p||2(^+l-m) */m2mn 

^ 2m!/ + 1 - m!2^+i-™ " " ^ ' 

m=0 
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□ 

Finally, the next proposition is used to evaluate a certain sum of normally ordered monomials in 
a{Oc'^t) and a*{Ocft), appearing in the computation of high moments of the random variable Ot- 

Proposition 8.3. For any £ gN, F € L^(M^), Vi)V'2 G J^, we have 
e 



m=0 



Proof. The proof proceeds by induction over ^ E N. For i = 0, the claim is clearly correct. We 
assume now it holds for i — 1, and we prove it for i. To this end, we use that, for every k > 1, 

ia{F) + a*iF)) : (a(F) + a*(F))'= :=: (a(F) + a*(F))^+i : +k\\Ff : (a(F) + a*iF))'''^ : (8.5) 

Applying the last identity twice, we conclude that 

(Vi, : (a(F) + a*{F) f : ^2> = (Vi, {a{F) + a*(F))2 : (a(F) + a*{F)f^'-'^ : ^2) 

- (M - 3)(V'i, : (a(F) + a*{F)f<^'-'^ : ^2) (8.6) 

- (2^ - 2)(2^ - 3)(V'i, : (a(F) + a*(F))2(^-2) : V2) 

Using the induction assumption, we find 

(V'l, (a(F) + a*(F))2 : (a(F) + a*(F))2(^-i) : ^^2) 
= (Vi,(a(F) + a*(F))2V2) 

m=0 

Again, using (jS.Sp . we obtain 

(V^i, (a(F) + a*(F))2 : (a(F) + a*(F))2{^-i) : V^a) 
= {i^i,{a{F) + a*{F)f%) 

m=0 
m=0 

e-2 



-4 2(^-l)!||Ff(^+^-) _ „VF112(-1) • \ 

2^^ 2(,„-l)!^-l-,„!2^-i-™^^i' • + " • 
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Combining the last equation with (|8.6p , we deduce that 

(V'i,:(a(F) + a*(F))2^:^2> 

= {i;,,{a{F) + a*{F)f'i^2) 

m=0 
m=0 

- 1^ 2(m-l)!£-l-m!2^-i-^^^' " ^"^^^ + " ^^^^ ' 

m=l 

Changing summation variables, we find 
= (V'i,(a(F) + a*(F))2V2) 

m=l 

_ g 2(.-l)!(4,„ + l)||Fr-) „ 

771=0 

m=0 

= (V^i,(a(F) + a*(F))2V2) 

_ ((£ _ i)(2^ _ 3) + 4£ _ 3) ||i7||2(^^^ : + : V'a) 

2(^-1)! 2(£-l)! 
(£-1)12^-1 ^ (^-2)!2^-2 

m=l 

X (2m(2m - 1) + 2(4m + 1)(£ - m) + - m){l - m - 1)) 



Since 



2/' 

l)(2£-3)+4^-3 = ^(2^-l) 



(2(£-l)!2' 



2(£-l)! ^ 2(^-1)! _ 2i\ 



(£-l)!2^-i (£-2)!2^- 



and 



2m(2m - 1) + 2(4m + 1)(£ - m) + 4(£ - m)(£ - m - 1) = 2^(2^ - 1) 
we conclude that 

(Vi, : (a(F) + a*(F))2^ : ^2> = (V-i, (a(i^) + a*(i^))'V2> 



m=0 
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which imphes 



(^1, (a(F) + a*{F))'%) = , J (V^i, : HF) + a*(F))^™ : ^^2) 

m=0 

and concludes the proof of the proposition. □ 
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